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Abstract: We show how to simulate the equatorial section of the Schwarzschild metric through a flowing liquid crystal 

in its nematic phase. Inside a liquid crystal in the nematic phase, a traveling light ray feels an effective 
metric, whose properties are linked to perpendicular and parallel refractive indexes, rio e rie respectively, 
of the rod-like molecule of the liquid crystal. As these indexes depend on the scalar order parameter of the 
liquid crystal, the Beris-Edwards hydrodynamic theory is used to connect the order parameter with the 
velocity of a liquid crystal flow at each point. This way we calculate a radial velocity profile that simulates 
the equatorial section of the Schwarzschild metric, in the region outside of Schwarzschild's radius, in the 
nematic phase of the liquid crystal. In our model, the higher flow velocity can be on the order of some 
meters per second. 
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1. Introduction 

Research in analogue models explores the similarities among areas of physics as different among themselves as 
general relativity and condensed matter [1]. It maps the main characteristics from a subject to another bringing 
in phenomenological insights even if there is not an obvious connection between them. For example, defects 
in crystals can be described by an effective Riemann-Cartan spacetime [2, 3]. In fluid systems, Unruh's effect 
[4] shows the similarity between the phonon creation in a radial flowing fluid and Hawking's effect around a 
black hole. Also, many analogue models for different black hole's aspects have appeared in the literature [5]. 
These procedures admit, not only new understanding on the subjects involved, but the possibility to make well 
controlled experiments with systems , like dielectric fluids or liquid crystals, as well. In this paper we are 
particularly interested in these latter ones. 

Liquid crystals are materials composed by anisotropic molecules, like rods, having different phases between the 
liquid and crystalline one. In one of these, the nematic phase, the rod-like molecules on the average point towards 
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a specific direction, represented by the vector n [6, 7]. As shown by Kline and Kay [8], Fermat's principle can 
be used to calculate light's trajectories in inhomogeneous anisotropic media like liquid-crystalline materials. This 
says that, among all the possible paths linking the points A and B, the light will choose the one that minimizes 
the integral N{r)dl, where N{r) is the refractive index of light at the position indicated by the vector f and I 
is a parameter with the dimension of length. This computation of light's path can be substituted by a geometric 
approach [9] and it is equivalent to say that the light, within a liquid crystal in the nematic phase, does not feel 
the geometric properties of the laboratory, but feels the properties of an effective space, of an effective metric. In 
this geometric point of view, the light trajectories are geodesies in the effective geometry. This makes a connection 
between nematic liquid crystals and space-time. For example, hedgehog-like defects in the nematic phase simulate 
the global monopole space-time [9], while the disclination-like defects simulate the metric of a cosmic string [9, 10]. 
The geometric properties of the "space-times" corresponding to the defects in liquid crystals depend on the 
perpendicular and parallel refractive indexes of the liquid crystal molecule in use, Uo e respectively. As these 
indexes are temperature dependent, it's possible to describe them in terms of a scalar function that measures 
how ordered is the nematic phase: the scalar order parameter q. If, for different points of the liquid crystalline 
material, there are different values of local velocity, then the Beris-Edwards hydrodynamic theory [11] can relate 
the local velocity of the liquid crystalline bulk to the order parameter q. This way, different velocity profiles 
generate different refractive index dependent functions and, thus, new effective metrics. 

With these ideas in mind, we study the case where light is propagating in the plane z = cte. of a configuration 
like a (fc = l,c = 0) disclination defect. This defect is a cylindrically symmetric linear one generated as follows. 
We assume that the liquid crystal flows radially inwards on a flat surface with a drain as in a sink. This tend to 
align the molecules radially in the conflguration known as radial disclination with unit topological charge [7, 9]. 
Far from this conflguration, the disordered phase is recovered. When a velocity fleld is applied in the bulk, the 
molecules orient according to the velocity gradient, generating an orientational order and a non-null value of the 
scalar order parameter q, that measures the strength of the ordering. Therefore, we calculate the velocity proflle 
v{r), that is the velocity of the liquid crystalline bulk at the distance r from the core of the conflguration. By 
doing this we simulate the equatorial section, 6 — tt/2, of the Schwarzschild metric: 

ds^ = - (^1 - dt^ + -jyKm^ + + sill' ^d<l,') . (1) 

This work is presented in the following way. First, we show the effective metric that describes the light propagation 
around a (A; = l,c = 0) disclination defect and how this propagation can be described by the order parameter 
q of the liquid crystalline material. Afterwards, we consider the liquid crystal flowing radially and we use the 
Beris-Edwards theory to analyze the dependence of the order parameter of the material with the flowing velocity 
module. In these two cases we consider the more general situation of three space dimensions. Finally, we employ 
the result from the second part in the flrst and we compare with the Schwarzschild metric written in isotropic 
coordinates [12]. 
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2. Geometric Model 

The analogy between Fermat's principle and the calculus of null geodesies in a Riemannian manifold has been 
known for a long time [13]. This relationship can be represented [9] by the following equation when we treat the 
extraordinary light [14] paths : 

being N^{r) the refractive index of the extraordinary light. In the (fc = l,c = 0) disclination case, the effective 
metric is [9] 

ds^ = b'^df^ + f'^dcj)^ + dz^, (2) 

where h = f = n^r and remembering that Uo and rie are respectively the refractive indexes of the liquid 
crystal molecule. 

Theses indexes can be written as temperature functions by [15] 



n.iT)^A-BT-^-^{,-L 



(3) 



Here, A, B, /3, (An)^ and Tc are constants that characterize the material, where (An)^ is the birefringence of 
the liquid crystal in the crystalline phase and Tc is the transition temperature between the nematic and isotropic 
phases. Considering the thermal dependence of the order parameter q through Haller's approximation [16] 

.(r)^(i-- 

the equation (3) changes to 

n4.) = A-BT.(l-gV.)-i^ 
ne {q) = A- BTc (l - q^'"^ + 

and thus the equation (2) becomes 

'a - BTc (l~q 



3 ' 
2(An)og 



ds^ 



1/I3\ _ ^(An)o9 2 



df^ + f^ddP + dz^ . 



yl_BT,(l-gi//9) + H(^^ 

The theory developed in [9] does not take into account the time coordinate but this coordinate must be taken 
into account in further calculations. In fact, the time coordinate enters the equation above in the simplest way. 
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since the geodesic equations of a manifold with the metric ds^ — A{r)dr'^ + [B{r)d6^ + C!{r) sin^ OdcfP^ are the 
null geodesies of the manifold with the metric ds'^ = —dt^ + A(r)dr'^ + [B{r)d6^ + C(r) sin^ , resulting in 

2 



A - BT, 



3 



A - BTc (1 - gi/") + 



2(An)og 
3 



+ f'^dS'^ + dz^ . 



We will back to equation (4) later. 



df 



(4) 



3. The Beris-Edwards Theory 

The relationship between the scalar order parameter g and the velocity v{f) of the liquid crystalline bulk on the 
position r is expressed by the Beris-Edwards equation [11] 



(9t+i;- V)Q-S(W,Q) =rH, 



(5) 



where Q is the order parameter tensor and F is a collective rotational diffusion constant. The tensor S (W,Q) 
depends on the gradient velocity tensor Wafi — dpVa as 



S(W,Q) = (^D + !^) (Q + 1/3) + (Q + 1/3) (CD - fi) 
- 2C(Q + I/3)Tr{QW}, 



(6) 



where ^ is the aspect ratio of the liquid crystal molecule, D — (W-|-W"'')/2 and Q — (W — W^)/2, the symmetric 
and antisymmetric parts, respectively, of the tensor W. The molecular field H on the right-hand side of equation 
(5) depends on the derivative of the free energy T by [17, 18] 

'ST] 



H 



-g + (I/3)T. 



5Q 



(7) 



= aQ + 6 (Q - (I/3)Tr [Q^] ) - hTr [Q^] 
+ kV^Q. 

In this section, we apply the Beris-Edwards equation (5) in the case where the liquid crystal bulk is moving 
radially towards a cylinder with radius ro 7^ (in microscopic scale) and the initial director is n = f , in cylindrical 
coordinates. This setup has the purpose to mimic the director of the (fc = l,c = 0) disclination, whose the 
effective metric we already know. Moreover, as the bulk flows continually inward, it is important that at r = ro 
there are sinks that removes the liquid crystal from the bulk, avoiding that the continuity law be violated. Thus, 
we use n = f to describe the order parameter tensor [Qij — Uirij — ^Sij) with the velocity dependence 

/ 2 n n \ 



Q{v{r}) = q{v{r)) 





-i 

-i 



(8) 



3 / 
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where v — \v\ = \v{r)\. As shown m [6], q{v) oc |Vt;(r)| and the cyhndrical symmetry of our setup says that the 
velocity v doesn't depend on the (j> and z coordinates, allowing one to write q{v{f)) = q(r). Therefore the flow's 
direction determines just where the maximum and the minimum of q occur. Expressing v{r) — v{r)r, the tensor 
S (W,Q) becomes 



S (W,Q) = 2^q (r) dr [v{r)] 1 - 



Mr) 



S ^ 


^ ^ 



Recalling that q{v) = q{r) and substituting (8) in (7), we obtain the molecular field H 



(9) 



Han = H,, = 



2aq{r) 2q\r)h 4:hq''{r) 
3 ^ 9 9 
2Kr^^dr [rdr [q{r)]] 



+ 



aq{r) q^{r)b 2hq-\r) 
~ 9~ ^ 9 

nr^'^dr [rdr [q{r)]] 



(10) 



The remaining terms are null. 

Substituting (8), (9) and (10) in (5), we get a matrix equation, whose solution rules the relationship between the 
scalar order parameter and the modulus of the velocity 



q{r) = 



2 + 



3C ' 



(11) 



where C is a constant defined by boundary conditions. 

4. The Effective Schwarzschild Metric 

When a metric is described by isotropic coordinates, the spatial section of the new metric is flat and all of it is 
multiplied by a function. Following the steps presented in [12], we rewrite a metric similar to (4) in the plane 
z = const., 



ds^ — —dt^ + p'^{f)dr^ + f'^dtf)^ , 



(12) 



in terms of isotropic coordinates. So (12) becomes 



ds^ = —dt^ + (p) [d/f + p^d(fP^ , 



(13) 
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where A (p) is a function to be determined and p is a new radial coordinate. Comparing tlie radial and angular 
terms from the last two equations we get 

p(f)df = \{p)dp, 

(14) 

f = A (p) p. 

For the metric (13) to have null geodesies that coincide with the null geodesies of the Schwarzschild one (1), i.e., 
for the metric (13) to be conformally related or conformal to (1) , we find in [5, 12] that 



1 + ^ 



3 



where M, in terms of the liquid crystal, represents the strength of the cylinder, that creates the disclination- 
like director, to divert the light. Substituting (15) in (14) and considering h'P <g 1, since we are dealing with 
microscopic ro, we obtain expressions for p{f) and p that allow (12) to have the same geodesies as the null 
geodesies of the Schwarzschild metric. We get 



where 



ds^ = -dr + ( 1 - —j dr + rd(j), (16) 



54f^M + 2f-3Af — - — 



fM 



8f3 - 54f2Af 



The last result has the following interpretation. For the Schwarzschild metric to be represented by the equation 

(4), the square root of the function associated with the radial coordinate on the equations (16) and (4) must be 

1 

equal. Doing this equality and noticing that ^ 1 (for ordinary values of /3 [15]), we find the radial dependence 
of q to be 



q{f) = 



48Mf' {A + BTc 



(An)o 
82Mf^ +2f= - 9M]"^ , 



X [32f ^ - 432Aff* + 8f ^ 

(17) 



where again we considered Al'^ ^ 1. This is the necessary scalar order parameter that must exist in a (fc = l,c — 0) 
disclination-like configuration in order to the metric (4) has geodesies that coincides with the Schwarzschild null 
geodesies for the consideration A/^ <C 1, that is equivalent to the regime of weak-field approximation. Notice that, 
up to now, we have obtained the most general expression for the scalar order parameter q, since no consideration 
was made on what is creating this specific order parameter (flowing bulk, temperature and/or electric field 
gradient, etc.). Now, considering that we are dealing with a fiowing liquid crystal, described by the case in the 
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previous section, and substituting equation (17) into (11), we find tlie velocity profile that creates the order 
parameter (17) to be 

w(f) = { [U4CMf^ {A + BTc)] X [(6f^ - 1296Mf* 

+24f ^ - 246A/f ^ + 6f - 27A/) (An)^ (18) 
- 96Mf^ {A + BTc) }^ . 

However, when this bulk's velocity is close to fight's velocity, drag effects must be considered [5]. 
It is important to note that we could use the Gordon metric [5] to determine a velocity profile that simulates 
Schwarzschild's metric, since we are treating with a fiowing medium and using the refractive index to calculate the 
behavior of light. However, the Gordon metric is used only when the refractive index of the medium is constant 
and here we are dealing with a varying refractive index, since [7, 9] 

Neir) = rio cos^ q + sin^ a, (19) 

where a is the angle between the propagation of light and the director n, and Uo and Tie are no longer constant, 
because they depend on the scalar order parameter q. 

Summarizing, every time we have a liquid crystal fiowing toward a disclination, as described in the previous 
section, and the bulk velocity is given by (18), any light ray traveling in this background will behave as in the 
presence of the Schwarzschild metric (1). But the key point is the equation (17), that shows how the order 
parameter must be, on the conditions above cited, for the Schwarzschild metric to be simulated. 

4.1. An Example 

Let's use the informations about the 5CB liquid crystal [15, 19] and apply them in equation (18) to obtain the 
velocity profile that simulates the Schwarzschild geometry. 

We consider that the liquid crystal system has a disclination-like configuration with a finite size core, but yet 
microscopic, and the liquid crystal fiows inward this core. Using the possible values of order parameter when the 
temperature of the liquid crystal varies on the nematic phase [7] , we admit that the maximum value of the order 
parameter q{r) is 0.5. Remembering this, we feed the equation (17) with the 5CB data and, admitting that the 
maximum value for the order parameter is reached on the surface of the defect (in this example, the core measures 
ro = 3 ■ 10"^ m ^ f ~ 5.12 • 10"^), the equation (17) provides us M ^ 1.13 • 10"^. 

It's important to mention that M is related to the Schwarzschild's radius, fs = 2M ~ 2.26 • 10"^ m < ro = 
3 • 10^''' m, and how the light travels outside this radius, following an open path. Also, we notice that increasing 
ro or the value of the order parameter on the core's surface, makes M increase. 

Having M on hands, the substitution of (17) in (11) gives C. However, this constant has an arbitrariness on 
its determination, because C depends on the boundary conditions. Admitting that the bulk's velocity as the 
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Figure 1. Velocity profile that simulates the Schwarzschild's metric around a disclination-like configuration with radius tq — 
3 ■ 10"^ m. 



fluid reaches the drain to be 1 m/s, we get C ~ 0.334. Finally, using all the constants found up to now in (18), 
we generate FIG. 1 which represents the velocity profile that simulates Schwarzschild's metric. Because of the 
arbitrariness of C, any other value of this constant just shifts up or downward the points in this graph, holding 
the behavior of the velocity's gradients. 

Now, let's evaluate the value of the fluid's velocity on the cylinder's surface that simulates the Schwarzschild 
metric obtained by Gordon metric. The algebraic expression for the flow velocity at a distance r is [20] 

VGordon{r) = —==^===, (20) 

where n is the refractive index of the medium. The latter quantity should be A^e since we are dealing with 
extraordinary light ray. Thus, equating n to (19) for 5CB on the cylinder's surface and for an inward light ray 
(a = 0), we have n = 1.54. Substituting this value, AI — 1.13 • 10^'' and r — ro = 3 ■ 10^'' m on (20), we obtain 
VGordon = 1.79 • 10* m/s, a value 10* times higher than one we obtained with C = 0.334, v — 1 m/s. Therefore, 
the use of the Gordon metric in this problem is experimentally infeasible. 

It's worth to note that the light in this flowing system describes open paths. Following the steps to calculate the 
light path as shown in most standard General Relativity books [12, 21], the values of ro and A4 generate open 
trajectories. Thus, we are not concerned with investigating the existence of a horizon. 
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5. Conclusion 

We have presented how a hght beam in a special configuration of a flowing nematic hquid crystal can effectively 
behave as if in the presence of a Schwarzschild space-time. Using the equivalence of Fermat's principle and the 
existence of an effective metric [9], Haller's approximation allows us to use the Beris-Ed wards' theory [11] to relate 
the properties of the effective metric to the local velocity of the liquid crystalline bulk. After that, Schwarzschild's 
metric described in terms of isotropic coordinates [5, 12] permits us to use a cylindrically symmetric configuration 
of the nematic liquid crystal and to obtain the appropriate radial velocity profile around this configuration that 
gives us Schwarzschild's metric. 

Remembering that we are dealing with velocities much lower than that of light, the distinguishing mark is the 
behavior of the bulk's velocity and not its absolute values. Thus, we could set velocity more appropriately within 
experimental possibilities. 

As a complement of this research, the study of interactive forces between the configurations here analyzed without 
the restriction of a total topological charge (otherwise, it's known as the dipole-dipole attractive force [22]) can 
be realized, as well as the theoretical study of some classical tests of general relativity [12] applied to this fiowing 
system. 

This work was partially supported by CNPq, CAPES, REDE NANOBIOTEC BRASIL and INCT-FCx. 
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